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We introduce a non-linear variant of the voter model, the q-voter model, in which q neighbors (with 
possible repetition) are consulted for a voter to change opinion. If the q neighbors agree, the voter 
takes their opinion; if they do not have an unanimous opinion, still a voter can flip its state with 
probability e. We solve the model on a fully connected network (i.e. in mean-field) and compute 
the exit probability as well as the average time to reach consensus by employing the backwards 
Fokker-Planck formalism and scaling arguments. We analyze the results in the perspective of a 
recently proposed Langevin equation aimed at describing generic phase transitions in systems with 
two (Z2 symmetric) absorbing states. In particular, by deriving explicitly the coefficients of such a 
Langevin equation as a function of the microscopic flipping probabilities, we find that in mean-field 
the g-voter model exhibits a disordered phase for high e and an ordered one for low e with three 
possible ways to go from one to the other: (i) a unique (generalized voter-like) transition, (ii) a series 
of two consecutive transitions, one (Ising-like) in which the Z2 symmetry is broken and a separate 
one (in the directed percolation class) in which the system falls into an absorbing state, and (iii) 
a series of two transitions, including an intermediate regime in which the final state depends on 
initial conditions. This third (so far unexplored) scenario, in which a new type of ordering dynamics 
emerges, is rationalized and found to be specific of mean-field, i.e. fluctuations are explicitly shown 
to wash it out in spatially extended systems. 

PACS numbers: 05.70.Ln, 02.50.2r, 05.50.+q, 64.60. Ht 



I. INTRODUCTION 

In a situation where one has to choose between two 
alternatives that appear equally agreeable, a rather com- 
mon way to remove the uncertainty is to copy what some- 
body else (randomly selected among acquaintances) does. 
The voter model dynamics [l[ describes precisely this sit- 
uation: Agents placed on the vertices of a graph are char- 
acterized by a binary (spin) variable ±1; at each time 
step two nearest neighbor vertices are selected and the 
first copies the state of the second. This can cquivalcntly 
be expressed by the flipping probability f(x): A vertex 
with a fraction x of disagreeing neighbors has a linear 
probability f(x) = x to flip. The iteration of this sim- 
ple (parameter free) rule gives rise to nontrivial ordering 
phenomena, that have drawn the attention of many schol- 
ars, both in physics Q and beyond 0]. From the point of 
view of statistical physics, voter dynamics stands out as 
one of the very few nonequilibrium processes amenable 
of exact analytical treatment in any dimension Q. In 
more physical terms, it owes its special character to the 
absence of surface tension [j| . Contrary to the more com- 
mon curvature-driven dynamics [||, in voter dynamics 
curved interfaces do not tend to reduce their curvature 
and assume a straight shape. This induces a slow do- 
main growth characterized, in two dimensions, by a log- 
arithmic decay of the density of active links (i.e. links 



connecting sites with opposite opinion states). 

A natural and relevant question, that has attracted in- 
terest in the past years, has to do with how generic the 
voter behavior is. Early work already pointed out 

that small changes in the dynamics destroy the voter be- 
havior in two dimensions. The voter model turns out to 
sit at the transition between a ferromagnetic and a para- 
magnetic phase but it is only a point in a generalized 
parameter space and any perturbation leads to a drasti- 
cally different behavior. These and other results poin ting 
towards the fragility of the voter behavior [lCJ, [HI, LL2|, LL3[ 
raise the question of whether the voter model is a pe- 
culiar exception or the representative of a more generic 
class of models. 

An answer to this question was provided by Dornic et 
al. Q. The authors of this work showed explicitly the 
existence of models, different from the pure linear voter, 
that nevertheless exhibit its typical dynamical features. 
This led them to conjecture that there is a proper general- 
ized voter (GV) universality class encompassing systems 
at "an order-disorder transition driven by the interfacial 
noise between two absorbing states possessing equivalent 
dynamical roles, this symmetry being enforced either by 
Z2 symmetry of the local rules, or by the global conser- 
vation of the magnetization" Q. 

Further progress in the understanding of this issue has 
been made by Al Hammal et al. [Til ], who worked out 
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a generic Langevin equation [15j for critical phenomena 
with two symmetric absorbing states and identified con- 
ditions for having a transition from order to disorder be- 
longing to the GV class. Note that in voter- like models 
there are two different competing phenomena: One is the 
breaking of the Z2 symmetry and the other one is the pos- 
sibility for the system to get trapped into an absorbing 
state. If both occur in unison then the transition point 
is in the GV class. Instead, if they occur separately, the 
Z2 symmetry is broken first (i.e. an Ising like transition 
occurs, and the system changes from paramagnetic to 
ferromagnetic) and afterwards the system falls into the 
corresponding absorbing state (i.e a directed percolation 
like transition) In this sense, the GV class can be 

rationalized as the superposition of Ising and directed 
percolation phase transitions. The fact that the voter 
transition can be split into two different ones was first 
reported in [T^ . 

The picture devised in Ref. [3| on the basis of generic 
symmetry arguments has been recently substantiated by 
Vazquez and Lopez fl7l |. Starting from the microscopic 
spin dynamics of a non-linear voter model, they have de- 
rived an explicit Langevin equation for the magnetiza- 
tion, that coincides with the one conjectured in Ref. [Til ]. 
In this way, it is possible to precisely determine, depend- 
ing on the analytical form of the microscopic flipping 
probability f(x), which of the two scenarios above oc- 
curs. 

In this paper we provide an assessment of the picture 
presented in Refs. \l4 . [l7j by proposing a microscopi- 
cally motivated non-linear voter model and analyzing it, 
both at the mean-field level and numerically. The model 
we consider, the (/-voter model, represents a simple gen- 
eralization of the original voter model: each individual 
interacts with a set of q of his nearest neighbors; if all q 
neighbors share the same state, the individual conforms 
to this state. Otherwise, if the q neighbors do not agree, 
the individual flips with a probability e. The (/-voter 
model is directly inspired in models of ordering dynamics 
in which each update step involves more than two indi- 
viduals [H, [SEIUS [UIUlj an d represents a practical 
and simple way to introduce nonlinearity in the voter 
dynamics at a microscopic level. 

We study the model phenomenology analytically via 
a mean-field approach, and numerically in two dimen- 
sions, uncovering a rich phenomenology. Depending on 
the value of q, the model exhibits all the possible transi- 
tions of a system with two symmetric absorbing states, 
as described above. Interestingly, at mean-field level the 
voter behavior is restricted only to very specific values 
of q (q = 2 and q = 3), two separated phase transitions 
occur in 2 < q < 3, and, otherwise (</ < 2 and q > 3) we 
find apparently novel phenomenology (i.e. dependence 
on the initial conditions and a double transition of a dif- 
ferent type). Direct numerical simulations of the model 
on a fully connected network are presented to back up 
the mean-field results. On the other hand, in a d = 2 lat- 
tice, we recover the picture presented in Ref. [l4| . with a 



single voter- like critical point. 

The paper is structured as follows: After Section [HI 
where the model is defined at a microscopic level, in 
Sec. IIIII we perform a mean-field analysis, by means of 
both a Fokkcr-Planck and a Langevin approach; a more 
detailed study of the case q = 4 is also presented. Results 
for finite dimensional systems are presented in Sec. IIVI 
Finally, in Sec. [V] we discuss our findings. 



II. DEFINITION OF THE q- VOTER MODEL 

We consider a non-linear voter model defined on a lat- 
tice (or network) of N sites. Each site hosts a spin, with 
value ±1. The dynamics is given by the following update 
rule: 

• At a given time i, choose one spin at random, lo- 
cated at site i. 

• Choose at random q neighbors of site i. In order 
to simplify the numerical analysis, and allow for 
an arbitrary value of q in regular lattices, we con- 
sider here the possibility of repetition, i.e. a given 
neighbor can be selected more than once. 

• If all the q neighbors are in the same state, the 
original spin takes the value of the q neighbors. 

• Otherwise, if the q neighbors arc different, the orig- 
inal spin flips with probability e. 

• Time is updated t — ► t + 1/N. 

It is easy to see that this model is non-linear. Consider 
the probability that a site flips as a function of the frac- 
tion x of disagreeing neighbors, that is, 

f(x, q) = x q + e [1 - x q - (1 - x) q ] . (1) 

Notice that, although in the original definition of the 
model q is an integer, Eq. |1} makes sense for any q > 
and can therefore be considered as the definition of the 
q- voter model for real values of q. For 3=1, we recover, 
for any value of e, the standard voter model, namely 
f(x, 1) = x. Non-linear behavior arises for q 7^ 1. Ob- 
serve also that if the fraction of disagreeing neighbors 
vanishes, i.e. x = 0, the configuration is absorbing, 
/(0,?) = 0. 

The (/-voter model bears some resemblance to other 
opinion dynamics models introduced recently in the lit- 
erature. For example, the "vacillating voter" model [l8[ 
is very similar to the (/-voter with q = 2 and e = 1 , apart 
from the possibility to select twice the same neighbor 
(repetition). The case with q = 2 and e = is instead 
similar to the Sznajd model [19(, in its formulation where 
a pair of agreeing agents convince only one of their neigh- 
bors H. 
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III. MEAN-FIELD THEORY 



and a diffusion coefficient 



In order to gain an understanding of the g-voter model 
behavior, it is useful to consider it first at the mean-field 
level (that is, on a fully connected network), for which 
several analytical tools have been recently developed. 



A. Backward Fokker-Planck approach 

Following Refs. [13, HH, HH, we can study the mean- 
field theory of the g-voter model by applying the back- 
ward Fokker-Planck (BFP) technique Consider a 
time t, in which there are n spins in state +1. The state of 
the system is fully defined by this quantity, plus the tran- 
sition rates to go to a state with n± 1 spins in state +1. 
Denoting this transition probabilities by p n ±i. n , then 



Pn+l,n = (l-x)f(x,q) 
Pn-l,n = Xf(l - X, q) 

Pn,n — 1 Pn+l.n Pn— l,nj 



(2) 
(3) 
(4) 



all the rest of values of p n ',n being equal to zero. Here 
we consider x = n/N as the probability of selecting a +1 
spin when a vertex in randomly chosen, a simplification 
which provides valid results in the limit of large N. With 
this definition there are two absorbing states, n = and 
n = N (i.e. x = and x = 1). 

The quantity n performs in time a biased one dimen- 
sional random walk, between two absorbing states. The 
random walk is fully defined in terms of a backward mas- 
ter equation, taking the form (l5j 

dP(n,t\n',t') 
dt' 

= T(n' + l\n')[P{n,t\n' + l,t') - P(n,t\n' ,t')] 
+ T(n' -l\n')[P(n,t\n' -l,t')- P(n,t\n',t')], (5) 

where P{n, t |n', t') is the probability of having n spins +1 
at time t, provided there were n' at time t' < t. Eq. ([5|) 
is given in terms of the transition rates (per unit time) 



T(n\ri) 



Pn.n' 

~A~ ; 



(6) 



with A = 1/iV. The master equation can be transformed, 
via a diffusion approximation, into a BFP equation for 
the reduced variable x = n/N, by expanding Eq. (|5|) 
up to second order in A. In this expansion, the BFP 
equation takes the form 

dP(x,t\x',t') , ,,dP(x,t\x',t') 
- = v(x )- 



dt' 



dx' 

1 „ d 2 P(x,t\x',t') 
+ 2 D{X) d7 2 ' (7) 



with a drift 



v{x) 



A[T(n+ l|n) - 
(1 - x)f(x,q) 



T(n- 
xf(l 



l\n)} 

- x,q), 



D(x) = A 2 [T{n+ l|n) +T(n- l|n)] 

= jj[(l-x)f(x,q)+xf(l-x,q)). (9) 

For the generic BFP equation Eq. ([7|), the exit proba- 
bility E(x), i.e. the probability that, starting from an 
initial density a; of +1 vertices, the absorbing state +1 is 
reached, satisfies the differential equation [1 51 ] 



v(x)d x E(x) + -D{x)d 2 x E{x) = 0, 



(10) 



subject to the boundary conditions -E'(O) = and E(\) 
i . while the average time until consensus, T(N, x) is given 
by |H 



v{x)3 x T{N,x) + -D{x)d 2 x T(N,x) = -1, 



(11) 



with boundary conditions T(N, 0) = T(N, 1) = 0. 

The standard voter model in Eq. ([S]) corresponds to 
q = 1, for which we find v(x) = and D(x) = 
2x(l - x)/N. This leads to E(x) = x and T(N,x) = 
-A^[a;lna; + (1 - xY\n(l - x)]\ thus T(N, 1/2) - N and 
T(N, 1/N) ~ In N [24J ] . Moreover, it is easy to see from 
Eqs. (fTUl) and (fTT|) . that the condition v(x) = is nec- 
essary and sufficient to yield, for any diffusion D(x) (as 
long as D{x) cx ^), E(x) = x and T(N, 1/2) - N, that 
are the two main signatures of voter behavior in mean- 
field. 

In order to have mean-field voter behavior, we must 
then consider the cases in which the drift v(x) vanishes. 
Let's look at the different possibilities: 



• For q = 2, 



v(x) = (-1 + 2e)(l - x)x(l - 2x). 



(12) 



(8) 



Therefore, e = 1/2 leads to voter behavior. But 
for e = 1/2, f{x,2) = x, so that in this case the 
g-model coincides with the usual voter model. 

• For q = 3, 

v(x) = (-1 + 3e)(l -x)x{\ - 2x). (13) 

Again, e = 1/3 leads to zero drift and hence to 
voter behavior. However, in this case f(x, 3) = 
x 3 — x 2 + x ^ x, so that the 3- voter model is a case 
belonging nontrivially to the GV class. 

• For q = 4 instead, 

v(x) = (1- x)x(l - 2x) 

x [-l + 4e + a;(l-a:)(l-2e)]. (14) 

No value of e can cancel the drift; therefore, voter 
behavior is, in principle, not possible. 
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B. Langevin equation approach 

Further understanding is provided by applying the for- 
malism developed in Refs. [Til, [l7|. In this approach one 
focuses on the magnetization <f> — 2x—l. In this variable, 
the drift takes the form, at lowest level in powers of <j>, 



v(cj>) = (l-0 2 )(a0-60 3 ) 



(15) 



This corresponds to the usual terms in a continuous de- 
scription for systems with a Z2 symmetry (i.e in the Ising 
class [13); multiplied by a factor (1 — (f 1 ) imposing the 
existence of two absorbing states. The drift can be writ- 
ten as derived from a potential: v((f>) = —dV(4>)/d(j), i.e. 



a , 2 a + 6 4 



6 



(16) 



This function has 5 extrema, obtained from the condi- 
tion v{4>) = 0, which are <j> = 0, <fi = ±1 and <f> = ±v/f- 
Their role is clarified by the concavity of V ((/)), that turns 
out to be 



V"(0) 
V"(±l) 




-a, 

2(a-6), 



2a(l-- 



(17) 
(18) 

(19) 



The extrema at (origin) and ±1 (absorbing barriers) 
are always relevant. The extrema at make physi- 

cal sense only when < £ < 1, otherwise they are imag- 
inary or non-accessible. According to the interpretation 
in Ref. [l4[, there are the following possible scenarios, 
depending on 6 (see Fig. [T|): 

• For b > 0, if a < the system is paramagnetic; an 
Ising transition occurs for a = and (afterwards) 
an absorbing-state (directed percolation) transition 
takes place at a = b (see Fig. [TJt). 

• The case 6 = corresponds to the voter case in 
which the potential identically vanishes at the tran- 
sition point a = (see Fig. [Tp). 

• For b < 0, if a is very negative the system is para- 
magnetic; then, at a negative value, a = 6, a pair 
of symmetric new minima appear at ±1 (this gen- 
erates an "intermediate phase" with three compet- 
ing minima). At a = the stability of the origin 
changes (see Fig. [Tp), and only the minima at ±1 
remain. 

In Ref. [HI it was argued that the intermediate phase 
appearing for b < is absent in spatially extended sys- 
tems: fluctuations wash it away, and the central curve in 
Fig. [Tp becomes as the lowest one. The reason for this 
is simple: As soon as minima at the absorbing barriers 
appear, fluctuations become asymmetric, i.e. they can 
take the system from the minimum at the origin to the 
barriers, but not the other way around. Therefore, in 
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FIG. 1: (Color online) Potential V (</>), as defined by Eq. (TTSJ) 
for b > (a), b = (b), and b < (b). 



a "renormalized" picture the third case is argued to co- 
incide with the 6 = case, and thus lead to a unique 
GV transition. Hence, only two scenarios are expected 
to exist in the presence of fluctuations. 

The coefficients a and b can be explicitly computed 
for the q- voter by combining Eqs. (|15p . ([8]), and (fTJ for 
generic values of q and e. This leads to 



2- q+1 {q - 1) - 2e(l - 2- q+1 ) 



b = 2-«( 9 -l)(?-2)(l-! 
+ 2e[l-2-«(2- (7 + < ? 2 )]. 



(20) 
(21) 



Before discussing what occurs for generic values of q, we 
remark that the coefficient b can be made to vanish iden- 
tically only for q = 1, q = 2, and q = 3, while the 
coefficient a vanishes for any e only if q = 1. Hence, we 
recover the previous results in Sec. IIII At For q = 1 the 
potential vanishes; the g-voter model coincides with the 
usual voter (e plays no role and the system sits at a crit- 
ical point). For q = 2 one has 6 = and a = 1/2 — e: the 
potential is exactly zero for e = 1/2 so that one recovers 
voter behavior at the transition between a paramagnetic 
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FIG. 2: (Color online) Mean field phase diagram of the q- 
voter model. At q — 1 , q — 2 and q — 3 (marked with vertical 
lines) GV transitions occur. For any other value of q there are 
two different transitions: at si = (solid line) the up-down 
(Z2) symmetry is broken, while at £2 = (dashed line) the 
barriers at ±1 become absorbing states. For 2 < q < 3 there 
are an Ising transition followed by a directed percolation one. 
For q < 2 and q > 3 an intermediate phase exists. 



and a ferromagnetic phase. For q — 3, we have 6 = and 
a = 1/2(1 — 3e): again the potential vanishes for e = 1/3 
(transition point) and voter behavior is found, separating 
an ordered phase (for small e) from a disordered one. 

For analyzing what happens for generic values of q, 
it is useful to calculate the boundaries of the interval of 
values of e for which the extrema in ± have physical 
values. We define as £1 the value at which a = so that 
the extrema are in 0. From Eq. (|2H|) we obtain 



£1 



1 



2" 



(22) 



Instead, £2 is the value for which the stability at the 
origin changes, i.e. for which a = b. From Eq. (|2ip we 
find 



e 2 



2_ 

3 



2q 2 



2«+ 2 -2(4- q + q 2 ) 



(23) 



The behavior of E\ and £2 as a function of q is plotted in 
Fig. [2j It is important also to notice that, for q > 1, a is 
smaller than zero for large £ so that there is a paramag- 
netic phase above the black solid line and a ferromagnetic 
one for small e. This corresponds to intuition: e plays 
the role of a sort of noise in the dynamics. Instead, for 
q < 1 (not represented in Fig. \2§ the situation is reversed: 
the paramagnetic phase is now for small £ and the ferro- 
magnetic one is for large values of e. 

The nature of the transition between the two phases 
varies depending on the value of q. We have already 
established that q = 1 represents a voter line, while for 
q = 2, and q = 3 (marked with vertical lines in Fig. [SJ) 
there is a voter transition point for the appropriate values 
of £. 

For 2 < q < 3, £1 is larger than £ 2 , and the nature of 
the extrema is therefore as follows: 



• For £ > £1, is a minimum and ±1 are maxima. 
The model is then in the paramagnetic phase (case 
Al in Fig. [I]) . 

• For £2 < £ < £1, and ±1 are maxima, and 

are minima. The mode is in the ferromagnetic 
phase (case A2 in Fig. [T]) . 

• For £ < £2, is a maximum and ±1 are minima. 
The model is in the ferromagnetic absorbing phase 
(case A3 in Fig.[T]). 

Obviously, in this case the double transition scenario de- 
scribed in Ref. applies: Starting from large values of 
£, first (at £ = £1) a transition in the Ising class, from 
a paramagnet to a ferromagnet, occurs; then (at £1) a 
transition of directed percolation type appears and the 
system becomes fully ordered. The same scenario occurs 
for q < 1, where, as mentioned above, the paramagnetic 
phase is for e < E\ and the absorbing one is for e > £2. 

On the other hand, for 1 < q < 2 and q > 3, the 
relative positions of £1 and £2 arc swapped, £1 < £2. 
There is an intermediate interval £1 < e < £2 such that 
and ±1 are minima separated by maxima in ±yja/b. 
The nature of the extrema is then as follows: 

• For £ > £2, is a minimum and ±1 are maxima. 
The model is in the paramagnetic phase (case Bl 
in Fig. [U). 

• For £1 < £ < £2, and ±1 are minima, and ±'\/f 
are maxima (case B2 in Fig. [TJ. 

• For £ < £1, is a maximum and ±1 are minima. 
The model is then in the ferromagnetic absorbing 
phase (case B3 in Fig. [1]) . 

In the intermediate interval the (mean-field) system ex- 
hibits ferromagnetic or paramagnetic behavior depending 
on the initial condition, with basins of attraction deter- 
mined by the separatrices ± \J a/b. The transition is com- 
plicated and there is no voter behavior at a mean-field 
level, as will be illustrated in the forthcoming subsection. 



C. Analysis of the q — 4 case 

The case q = 4 is the smallest integer value of q for 
which the coefficients a and b cannot be made identically 
equal to zero simultaneously: 



Ue 



and 



1 



(24) 



Correspondingly we have £1 = 3/14 and £2 = 1/4. In 
Figure [3] we plot the average consensus time as a func- 
tion of N for q = 4 and several values of e, obtained by 
numerical simulations of the (/-voter model on fully con- 
nected networks of different size. For £ > 1/4 the growth 
is exponential, as expected in the paramagnetic phase. 
For £ < 3/14 the growth is logarithmic, as expected in 
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FIG. 3: (Color online) Average consensus time T(N,x — 1/2) 
as a function of TV, for a fully connected system with q — 4 
and different values of e. The straight line has slope 1/2. 
From bottom to top: e = 0.15, 3/14, 0.22, 0.23, 0.25, and 0.35. 



the ferromagnetic phase. For 3/14 < e < 1/4 there is 
a crossover, but asymptotically the growth is exponen- 
tial. For e = 3/14 the growth is proportional to N 1 / 2 , 
different from the voter linear behavior. 

A simple scaling argument allows us to understand the 
growth law N 1 ^ 2 for the average consensus time at the 
transition point e = 3/14. For q = 4 and e = 3/14 
the potential has parameters a = 0, b = — 1, hence 
V{4>) = — <7J 4 /4 at leading order. For initial conditions 
x = 1/2, = 0, so that at the beginning only diffusion 
matters. When <f> becomes sufficiently large drift comes 
in, the motion becomes ballistic and in an interval de- 
pending logarithmically on N consensus is reached. How 
much time is spent in the diffusion stage? This interval 
lasts a time t* that is estimated by equating the drift 
v ~ i' 
tion, 

During this time interval the diffusive motion gives a dis- 



6 3 with the e ffective velocity of the diffusion mo- 
V~Dt*/t* = s/D/t*. Hence t* = D/v 2 ~ D/(jP. 



placement <fi* = y/Dt* yielding t* = 
ing the two expressions gives D ~ (0" 
hand the diffusion coefficient D is D 
(1/N) 2 /(1/N) = 1/N, from which </>* 
the total time to consensus is 



\4>*) 2 /D. Equat- 
) 4 . On the other 
= (A$ 2 /(Ai) = 
- N- 1 / 4 . Hence 



T(N, 1/2) w ((f)*) 2 /D + \og{N) - N 1/2 + \og(N) - N 1/2 . 

(25) 

With respect to the exit probability E(x), we can apply 
the formalism in Sec. IIII Al solving Eq. (flTj|) for the drift 
and diffusion functions 



v(x) = ~-{l-x)x{l-2xf, 
1 



(26) 



D(x) = — (1 - x)x (24x 2 - 2Ax + U) , (27) 
to obtain, applying standard stochastic techniques p8j , 



E(x) 



fte-ZVy-W (24y 2 -24y+13)' 



dy 



/„ e -Tl(2y-i) 2 (2Ay 2 - 24y + 13) 



7^/72 



dy 



FIG. 4: (Color online) Exit probability E(x) as a function of 
x, and different system sizes N, for a fully connected system 
with q = 4 and e = 3/14: the larger the size the steeper the 
slope at x = 1/2. Full lines represent the numerical integra- 
tion of Eq. (|28[>; symbols stand for direct numerical simula- 
tions. 



In Fig. 0] we show the results of the numerical inte- 
gration of Eq. (|28|) for different values of N (solid lines), 
together with results form direct numerical simulations 
of the q- voter model in fully connected networks. It is 
apparent that E{x) tends to a step function at x = 1/2 
in the limit N — -> oo. No voter behavior is therefore 
observed in this mean-field case. 



IV. BEHAVIOR IN FINITE DIMENSIONS. 

A. d= 1 

Even if in d = 1 the number of different nearest neigh- 
bors (in a square-lattice) is 2, the parameter q can be 
kept arbitrary. It is easy to see that the one-dimensional 
q-voter model can be mapped onto the model of non- 
conservative voters recently introduced in Ref. [2{|. In 
such a model, the relevant parameter 7 is given by the 
ratio P2/P1, where pi is the flipping probability for a 
site surrounded by i disagreeing neighbors. In the q- 
voter model pi trivially equals 1 for any q, while p\ = 
2-i + e(l-2 1 -^). Hence 



2'/ 



1 + e(2<? - 2)' 



(29) 



(28) 



This equation leads to the conclusion that the value 
e = 1/2 yields, for any q, a voter behavior (7 = 2). 
Analogously one can see that no e can give the value 
7=1, implying the 'vacillating voter' behavior [l8| . The 
mapping from q and e to 7 allows us to deduce, from 
the results of Ref. 29;], the nontrivial shape of the exit 
probability E(x) in d = 1. 
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FIG. 7: (Color online) Numerical estimate of the average drift 
FIG. 5: (Color online) Plot of the inverse of the density of in the (/-voter model, as a function of the magnetization, (a) 
active links p(t) for different values of e. Parameters: q = 4, Square lattice, (b) Fully connected network. In this last plot, 
e = 1/4, L = 5000. the full line represents the expected theoretical value v(<j>) ~ 

(1- 



B. d = 2 

As in d = 1, here q is kept arbitrary. We have per- 
formed numerical simulations of the ordering dynamics 
of the q- voter model on a square lattice of size Lx L with 
L = 5000, for several values of q > 1. In all cases, we find 
a transition separating a paramagnetic phase for large e 
from a ferromagnetic one at low e. In order to investigate 
the nature of the transition we concentrate on the case 
q = 4, that in a fully connected graph yields non voter 
behavior. In Fig. [5], we plot the temporal behavior of the 
inverse of p, the fraction of active links in the system. 
At the critical point e = E\ = 1/4, 1/p grows logarith- 
mically, as expected for the voter universality class [1]. 
Analog results are found for other values of q (data not 
shown). Additional evidence proving that for e = 1/4 the 
q- voter model behaves exactly as the usual voter model 
is provided by measuring the correlation function C (r, t) . 
From the exact solution of the voter model in d = 2 [3(| 
it turns out that two different length scales are present 



4 _ 

S 

2- 





10-' 



in the system, leading to the nonstandard scaling form 

1 



C(r,t) = 



ln(16i) 



f[r 2 /(2t)], 



(30) 



where the f(x) = Ei(x) is the exponential integral func- 
tion [3l|. Figure [HI demonstrates that Eq. ([50)) is nicely 
obeyed by numerical simulations of the g-voter at the 
transition point for q — 4. This evidence, further con- 
firmed by the analysis of the exit probability, showing a 
linear behavior, leads to strong numerical confirmation 
that the scenario predicted by Al Hammal et al. [l4| is 
correct in d = 2: In finite dimensions, fluctuations renor- 
malize the deterministic potential so that the transition 
is in the GV class for any value of q. This rcnormalization 
effect can be directly observed by measuring numerically 
the drift term. From the Fokker-Planck equation, we can 
obtain the equation for the time evolution of the magne- 
tization, (0), namely [i~5| 



d{4>) 
dt 



(31) 



10 

r 2 /(2t) 



Therefore, a numerical evaluation of d{4>) j dt yields an 
estimate of (v(4>)). In Fig. [7] we plot the average drift as 
a function of the magnetization for the g-voter model in 
d = 2 and for the mean-field fully connected case. In the 
latter, (v((j))) shows a functional dependence compatible 
with the theoretical expectation v{<j>) ~ (1 — </> 2 )</> 3 for 
sufficiently large network size N . In the d = 2 case, on 
the other hand, fluctuations are able to quickly cancel 
the drift term, inducing thus an effective voter behavior 
in the limit of large N: b < renormalizes on large scales 
to b = 0, as predicted in (l4| . 



FIG. 6: (Color online) Collapse plot of the correlation func- 
tion C(r,t) for different times, according to Eq. (|30|l . Param- 
eters: q = 4, e = 1/4, L = 5000. 



V. CONCLUSIONS AND DISCUSSION 

We have introduced a non-linear variant of the voter 
model in which the opinion of q neighbors (with possible 
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repetition) is taken into account for a voter to change its 
own opinion. In particular, if all his q neighbors share 
the same state, an individual conforms to this state; oth- 
erwise, if the q neighbors do not agree, he flips with a 
probability e. Note that the model includes a noise effect 
controlled by e; still voters are not allowed to break the 
absorbing state condition and a consensus state remains 
indefinitely so. While the original definition of the model 
is meaningful only for integer values of q, analytical gen- 
eralization to arbitrary values of q, with q € [0, oo] is 
possible. In particular, after taking the continuous limit 
for the transition rates, q becomes a, not necessarily in- 
teger, parameter. 

We have studied the model analytically by applying 
a mean field analysis based on the backward Fokker- 
Plank formalism and the Langevin approach developed 
in Refs. [HI, [13] ■ These two approaches permit us to un- 
cover the rich and variate phenomenology of the g-voter 
model: 

• For q = 1 the model reduces to the standard voter 
model with exit probability proportional to x and 
average consensus time, starting from x = 1/2, 
growing linear with system size. 

• For q — 2, the model coincides with the voter one if 
e = 1/2. Instead, the system remains disordered for 
e > 1/2, or it orders exponentially fast for e < 1/2. 
Therefore the 2-voter model exhibits a "generalized 
voter transition". Note that our results also clar- 
ify the behavior of the Sznajd model. For q = 2 
and £ = the g-voter model practically coincides 
with Sznajd model, at least in the formulation of 
Ref. [HI , according to which one has to select a pair 
of neighbors and, if they are in the same state, an- 
other neighbor of the pair is set in the same state. 
From this point of view, Sznajd model is just a fer- 
romagnetic model in its ordered phase. 

• For q — 3 there is a voter like transition at e = 1/3, 
separating two phases as those described for q = 
2 but, contrarily to the cases before, the flipping 
probabilities are non-linear: the 3- voter model is an 
example belonging non-trivially to the generalized 
voter class. 

• For 2 < q < 3 there is no voter transition. Instead 
the system experiences a sequence of two transi- 
tions. Starting with large values of e and reducing 
it progressively, first the Zi symmetry is broken 
(Ising transition) and afterward the system orders 
into an absorbing state at a directed-percolation 
like transition. 



• For 1 < q < 2 and q > 3 the mean-field approach 
predicts a non- voter transition, characterized by an 
exit probability with a Heaviside O-function shape 
and a consensus time which increases with system 
size as N 1 / 2 rather than linearly. 
All these results have been verified in numerical simu- 
lations of the model on a fully connected lattice. On the 
other hand, in spatially extended systems this last (third) 
scenario does not appear, as predicted by Al Hammal 
et al. fLU ]. The reason for this is, as we have numeri- 
cally verified, that fluctuations wash out the intermedi- 
ate regime in which three stable states exist. Indeed, as 
fluctuations can take the system from the origin to any 
of the absorbing states but not the other way around, 
effectively, on sufficiently large scales, the stable state at 
the origin plays no role, and the system exhibits a single 
ordering transition in the generalized-voter class. 

Previous results were obtained under the rule that, 
among the q neighbors involved in the dynamical step, 
each given neighbor may be selected more than once. If 
the possibility of repetition is explictly forbidden, mean- 
field results clearly do not change, but in finite dimen- 
sions variations are possible. We have performed numer- 
ical simulations in d = 2, showing that for q < 4 no 
qualitative change occurs: there is a GV transition be- 
tween a disordered phase for large e and an ordered one 
for small e. Things change only in the special case q = 4, 
where all neighbors of a site are considered. The ordered 
phase disappears and for any value of £ a disordered state 
is reached. 

In summary, the g-voter model is a simple non-linear 
extension of the voter model exhibiting a rich and inter- 
esting phenomenology and illustrating how apparently 
innocuous changes in the microscopic dynamics can lead 
to different types of collective phenomena, and in partic- 
ular to different paths to reach consensus. 
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